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Abstract 

Back in the early days of polyomino enumeration, a model called 
column-convex polyominoes was introduced and its area generating func- 
tion was found. That generating function is rational: the numerator 
has degree four and the denominator has degree three. Let a column- 
duplex polyomino be a polyomino whose columns can have either one or 
two connected components. In spite of being an immediate generaliza- 
tion of column-convex polyominoes, column-duplex polyominoes do not 
seem amenable to any kind of exact enumeration. Hence, we are inter- 
ested in models which interpolate between column-convex polyominoes 
and column-duplex polyominoes. This paper is concerned with a model 
which we have called simplex-duplex polyominoes. A simplex-duplex poly- 
omino is a column-duplex polyomino in which there is no occurrence of two 
adjacent columns each having two connected components. Simplex-duplex 
polyominoes are not easy to deal with, but their area generating fuction 
can still be found. To find this generating function, we use an upgraded 
version of the Temperley method. Though that technique is widely used 
in these times, our application presents two interesting features. Firstly, 
we add one or two columns at a time, thus bypassing those simplex-duplex 
polyominoes which end with a two-component column. (It is somewhat 
more usual to add just one column at a time.) Secondly, we obtain a 
functional equation that involves both the first and the second derivatives 
of the sought-for generating function. (Such equations usually involve the 
first derivative only. In some cases, no derivative is involved at all.) Right 
because of this latter interesting feature, the Temperley method produces 
a very complicated formula for the generating function. Anyway, from 
that formula it is easy to compute Taylor polynomials. Thus we get 
plenty of evidence that the number of n-celled simplex-duplex polyomi- 
noes behaves asymptotically as 0.119443 • 3.522020™ . For comparison, the 
number of n-celled column-convex polyominoes behaves asymptotically as 
0.180916 • 3.205569™. 
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1 Introduction 



Enumeration of polyominoes is a well-established topic in mathematics, physics 
and chemistry It would be ideal to know how many polyominoes can be 
made up of n cells, or how many borders of polyominoes can be made up of 
m unit edges. Realistically speaking, however, neither of these questions will 
be answered in the foreseeable future. So it is a sound idea to confine atten- 
tion to sets of polyominoes which do not contain all polyominoes, but anyway 
contain rather many of them. Sets of that kind sometimes allow of the enu- 
meration by perimeter and/or by area. Yet, in defining such a "friendly" set of 
polyominoes, there is not much room left for poetic license. That kind of defi- 
nitions almost always require that the elements of the set be either in a certain 
sense convex or in a certain sense directed. A classical example are column- 
convex polyominoes: their definition requires that the intersection between the 
polyomino and any vertical straight line must be a connected set. The area 
generating function for column-convex polyominoes has been known at least 
since 1956 [131 114] . The perimeter generating function for column-convex poly- 
ominoes is also not a recent result [H [S] El- However, there exists an ongoing 
line of research aimed at constructing models which are supersets of column- 
convex polyominoes and are still solvable. The models constructed so far are 
called m-convex polygons (12] . prudent polygons [10] , cheesy polyominoes [7], 
polyominoes with cheesy blocks [8], and column- subconv ex polyominoes [9]. The 
former two models allow of the enumeration by perimeter, whereas the latter 
three models admit of the enumeration by area. (Thus, it is still an open problem 
to think out a model which significantly generalizes column-convex polyominoes 
and admits of both the enumeration by perimeter and the enumeration by area.) 

Let a column-duplex polyomino be a polyomino whose columns can have ei- 
ther one or two connected components. Somewhat surprisingly, column-duplex 
polyominoes are already too hard a model: they seem to admit of no exact enu- 
meration. This paper is concerned with a model which we have called simplex- 
duplex polyominoes. A cheesy polyomino, a polyomino with cheesy blocks, a 
column-subconvex polyomino and a simplex-duplex polyomino are all column- 
duplex polyominoes. In the first three models, the definition requires (perhaps 
among other things) that the gap within a two-component column must not 
be greater than m cells in size, where m is a positive integer which we fix in 
advance. In the fourth model (simplex-duplex polyominoes), the gap within a 
column can be of any size, but it is not allowed that two two-component columns 
come one immediately after the other. 

The main result of this paper is a formula for a bivariate generating func- 
tion G(q,w), in which the coefficient of q n w is the number of simplex-duplex 
polyominoes having n cells and k two-component columns. To compute the for- 
mula for G{q 1 w), we make use of an upgraded version of the Temperley method 
[21 [14]. Our computations have some interesting features, as is described in 
the abstract above. Although G(q, w) and G(q, 1) prove to be very compli- 
cated g-series, from the formulas for G{q, w) and G(q, 1) it is easy to find Taylor 
polynomials of any reasonable degree. Thus we get plenty of evidence that 
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the number of n-celled simplex-duplex polyominoes behaves asymptotically as 
0.119443 • 3.522020™. This means that simplex-duplex polyominoes are one of 
the largest sets of polyominoes for which the area generating function is known 
(cf. [3]). We already enumerated some fairly large sets of polyominoes 0IH1III], 
but those polyominoes had hexagonal cells, whereas simplex-duplex polyomi- 
noes have square cells. However, it is fair to say that the best currently known 
lower bound on Klarner's constant is 3.980137 [1], which is a significantly greater 
number than our growth constant 3.522020. 

This paper continues as follows. In Section 2, we state the necessary def- 
initions and conventions. In Section 3, we establish and solve a functional 
equation. The function satisfying that equation is a generating function not 
for all simplex-duplex polyominoes, but only for those simplex-duplex polyomi- 
noes which end with a one-component column. Once the functional equation 
is solved, it is however very easy to obtain the area and duplex columns gener- 
ating function for all simplex-duplex polyominoes. The formula for this latter 
generating function is stated in Section 4. In the same section, there is also 
a corollary, in which we state the area generating function for column-convex 
polyominoes. In Section 5, there is a bit of asymptotic analysis. In Section 6 
we conclude, outlining further work prompted by our results. 

2 Definitions and conventions 

In this paper, a cell is a unit square whose vertices have integer coordinates. 
A plane figure P is a polyomino if P is a union of finitely many cells and the 
interior of P is connected. 

Let P and Q be two polyominoes. We consider P and Q to be equal if and 
only if there exists a translation / such that f(P) = Q. 

Given a polyomino P, it is useful to partition the cells of P according to their 
horizontal projection. Each block of that partition is a column of P. With this 
definition, a column of a polyomino is not necessarily a connected set. However, 
it may happen that every column of a polyomino P is a connected set. In this 
case, P is called a column- convex polyomino. 

By a column-duplex polyomino, we mean a polyomino in which columns with 
three or more connected components are not allowed. Thus, each column of a 
column-duplex polyomino has cither one or two connected components. 

Let P be a column-duplex polyomino and let c be a column of P. If c has 
one connected component, we say that c is a simplex column. If c has two 
connected components, we say that c is a duplex column. 

A simplex- duplex polyomino is such a column-duplex polyomino in which 
consecutive duplex columns are not allowed. If c is a column of a simplex-duplex 
polyomino, and c is a (whether left or right) neighbour of a duplex column, then 
c must be a simplex column. See Figure 1. 

Let P be a simplex-duplex polyomino and let P have at least two simplex 
columns. Then we define the lower pivot cell of P to be the cell which is the 
right neighbour of the bottom cell of the second-last (i.e., second-rightmost) 
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Figure 1: A simplex-duplex polyomino. 
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among the simplex columns of P. We also define the upper pivot cell of P to be 
the right neighbour of the top cell of the second-last among the simplex columns 
of P. 

Let P be a simplex-duplex polyomino and let P have at least one duplex 
column. Then we define the lower inner pivot cell of P to be the right neighbour 
of the top cell of the lower component of the last among the duplex columns of 
P. We also define the upper inner pivot cell of P to be the right neighbour of 
the bottom cell of the upper component of the last among the duplex columns 
of P. 

Observe that the lower pivot cell of a simplex-duplex polyomino P is not 
necessarily contained in P. The same holds for the upper pivot cell, the lower 
inner pivot cell and the upper inner pivot cell of P. 

If a polyomino P is made up of n cells, we say that the area of P is n. 

Let S denote the set of those simplex-duplex polyominoes whose last (i.e., 
rightmost) column is a simplex column. It is convenient to first compute a 
generating function for the set S, and thence a generating fuction for all simplex- 
duplex polyominoes. So, let 

the height of the the number of 

A(q,t,w) = <f rea ° f P -t last c ° ,umn ° S P ■ w dupl " COIumns ° f P . 
Pes 

Next, we define three generating functions in two variables, q and w: Let 
A 1 = A(q,l,w),B 1 = ( ^=f=^ \ andd^-j ^y 1 ] } 

I ) with t—1 y ) with t—1 

For n G N, let T n be the set of simplex-duplex polyominoes having a duplex 
last column and an n-celled second-last column. Every element of T n is made 
up of three ingredients: 

1. the "left factor" consisting of all columns except the last column, 

2. the lower connected component of the last column, 

3. the upper connected component of the last column. 

The area and duplex columns generating function for ingredient No. 1 is 
(t n )A(q, t, w), where (t n ) stands for "the coefficient oft™". The area generating 
function for each of the other two ingredients is The three ingredients can 

be glued together in exactly ("g 1 ) = ("~ 1 K"~ 2 ) wa y Sj because there arc ("~ 1 ) 
ways to choose two non-adjacent cells in the last column of ingredient No. 1. 
Altogether, the area and duplex columns generating function for the set T n is 

[(t n )A(q, t, w)] ■ ("~ 1 K"~ 2 ) . (^j^j ■ w. As to all simplex-duplex polyominoes 
ending with a duplex column, their area and duplex columns generating function 
is 



E~ i W n )A{q,t, W )] ■ ( "' 1) 2 < "' 2) • (t^) 2 • " 
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By now it is clear that 



G = Ai + 




(1) 



(1-9) 



3 Establishing and solving a functional equation 

Henceforth the notation A(q, t, w) will be abbreviated as A(t). 

In order to obtain a functional equation for the generating function A(t), 
we are going to suitably partition the set S. (Recall that S is the set of those 
simplex-duplex polyominoes whose last column is a simplex column.) 

First, let S a be the set of those PeS which have no other simplex column 
than the last column. 

Let Sp be the set of those P £ S\S a which have the following two properties: 

• the second-last column is a simplex column, 

• the last column contains the lower pivot cell of P. 

Let Sj be the set of P G S \ S a having the following two properties: 

• the second-last column is a simplex column, 

• the last column does not contain the lower pivot cell of P. 

Thus, Sp U S 7 is the set of those P e S\ S a whose second-last column is a 
simplex column. 

Let Sg be the set of those P e S\S a which have the following three prop- 
erties: 

• the second-last column is a duplex column and the third-last column is 
(necessarily) a simplex column, 

• the lower component of the second-last column and the third-last column 
have no edge in common, 

• the lower pivot cell of P is contained in the upper component of the 
second-last column. 

Let S € be the set of P e S \ S a having the following three properties: 
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• the second-last column is a duplex column and the third-last column is a 
simplex column, 

• the lower component of the second-last column and the third-last column 
have no edge in common, 

• the lower pivot cell of P is contained in the hole of the second-last column. 

The definition of Sq is obtained from the definition of Ss by writing the 
word "upper" where the definition of Sg says "lower" , and by writing the word 
"lower" where the definition of Sg says "upper" . 

The definition of S v is obtained when the changes just described are made 
to the definition of S £ (instead of to the definition of Sg). However, since the 
set Srf appears in the proof of Theorem 1 below, it makes sense to define S. q 
explicitly: S n is the set of P e S \ S a having the following three properties: 

• the second-last column is a duplex column and the third-last column is a 
simplex column, 

• the upper component of the second-last column and the third-last column 
have no edge in common, 

• the upper pivot cell of P is contained in the hole of the second-last column. 

Thus, S s U S e U S{ U S v is the set of those P G S \ S a which, in ad- 
dition to having a duplex second-last column, also have the property that 
P \ (the last column of P) is not a polyomino. 

Let Sg be the set of P e S \ S a which have the following three properties: 

• the second-last column is a duplex column and the third-last column is a 
simplex column, 

• each of the two components of the second-last column has at least one 
edge in common with the third-last column, 

• the last column also has at least one edge in common with each of the two 
components of the second-last column. 

Let S L be the set of P £ S \ S a which have the following four properties: 

• the second-last column is a duplex column and the third-last column is a 
simplex column, 

• each of the two components of the second-last column has at least one 
edge in common with the third-last column, 

• the last column has at least one edge in common with the lower component 
of the second-last column, but docs not have any edges in common with 
the upper component of the second-last column, 
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• the last column does not contain the lower inner pivot cell of P. 



Let S K be the set of P £ S \ S a having the following four properties: 

• the second-last column is a duplex column and the third-last column is a 
simplex column, 

• each of the two components of the second-last column has at least one 
edge in common with the third-last column, 

• the last column has at least one edge in common with the lower component 
of the second-last column, but does not have any edges in common with 
the upper component of the second-last column, 

• the last column contains the lower inner pivot cell of P. 

The definition of S\ is obtained from the definition of S L by writing the 
word "upper" where the definition of S,, says "lower" , and by writing the word 
"lower" where the definition of S t says "upper" . 

The definition of S 1 ^ is obtained when the changes just described are made 
to the definition of S K (instead of to the definition of S L ). 

Thus, S e U S L U S K U S\ U 5^ is the set of those P £ S \ S a which, in 
addition to having a duplex second-last column, also have the property that 
P \ (the last column of P) is a polyomino. This means that S5 U S e U . . . U 
is the set of all P £ S \ S a whose second-last column is a duplex column. 

The sets S a , Sp, ... , form a partition of the set S. 

Let the notations A a (t), Ap(t), . . . , A^t) denote the parts of the series A(t) 
that come from the sets S a , Sp, . . . , S^, respectively. 

Theorem 1 We have: 



qt (1-5)2(1-^)3 ' 

Q t W 

Mt) = Mt) = (1 _ g) 3 (1 _ gt) 3 -^ 

q 5 t 3 w qH 2 w 
A € (t)=A v (t) - (l _ 3 )2(i _ qt y " Al ~ (1 - g)2(i - qt y • A W> ' 

Mt) = (l- g )»(l -*)■ •* 



(l-q)Hl-qty {\-q)\\-qt) 4 



Mt) = Mt) = {1 _^ 1 _ qty c 1 , 
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a m - a m qHw r qH3w r 



q^t^w qH w . , , 



(1-^)2(1-^)4 (l-g)2(l-gt)4 

Proof of Theorem 1. We shall only prove the formulas for A v (t), Ag(t) and 
A K (t). The other formulas are proved in a similar way. 
Consider this procedure: 

Procedure 1 Given a polyomino P £ S; nonempty columns C\, c 2 , c^, C5 and 
Cq, none of which has a hole; a column C4, which may be empty, but may not 
have a hole. First place the column c\ so that C\ lies two units to the right of 
the last column of P and so that the top of c\ is aligned with the top of the last 
column of P. Then place the columns c 2 , ■ ■ ■ , cq as shown in Figure 2. 

Procedure 1 does produce all elements of the set S v , but also produces some 
objects which do not lie in SL. Indeed, if the column ci is equally high or higher 
than the last column of P, we get an object in which the second-last and third- 
last columns do not have any edges in common. That object is obviously not a 
polyomino. To get rid of those undesired objects, we define a second procedure. 

Procedure 2 Given a polyomino P € S; nonempty columns c 2 , C3, C5 and cq, 
none of which has a hole; columns C\ and C4, which may be empty, but may not 
have a hole. First make a copy of the last column of P and place that copy so 
that the copy lies two units to the right of the last column of P and so that the 
top of the copy is aligned with the top of the last column of P. Then place the 
columns C\ , . . . , Cq as shown in Figure 3. 

What Procedure 2 produces is precisely the undesired part of the output of 
Procedure 1. Now, the generating function for the output of Procedure 1 is 



'l-qt (X-qf(X-qtr' Av (2) 

In the contribution of the polyomino P is A\, the contributions of C3 
and cq are each, the contributions of c\ , c 2 and C5 are is each, and the 
contribution of C4 is yz^j. The contributions are multiplied together because 
Procedure 1 specifies uniquely where each of the pieces is to be placed. 

The generating function for the output of Procedure 2 is 



In ©, the joint contribution of P and the copy of P's last column is A(qt), 
the contributions of C3 and cq are each, the contributions of c\ and C4 are 
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each, and the contributions of C2 and C5 are is each. Again, the 

contributions are multiplied together because Procedure 2 specifies uniquely 
where each of the pieces is to be placed. 
Subtracting ([3]) from ([2]) gives 

a , s q 5 t 3 w q A t 2 w 

A ^ = (1-qni-qt)* ■ Al (1-^2(1-^)4 ' 

which was to be demonstrated. 
Next, consider this procedure: 

Procedure 3 Given a polyomino P S S and five nonempty columns c\, . . . , C5, 
none of which has a hole. First choose a cell ( call it z) which lies in the last 
column of P and is not the top cell of that column. Then place the column c\ so 
that the top cell of C\ becomes the right neighbour of z. Then place the columns 
C2, . . . , C5 as shown in Figure 4- 

Procedure 3 produces all elements of the set Sg, but also produces some 
objects which do not lie in Sg. Fortunately, those undesired objects are easy to 
identify: they are precisely the elements of S^, the set whose generating function 
has just been found. Thus, we shall now compute the generating function for all 
objects produced by Procedure 3. Then we shall subtract A v {t) and the result 
will be A e (t). 

In Procedure 3, if the last column of P has n cells, the cell z can be chosen in 
n — 1 ways. Once z is chosen, the places where ci, . . . , C5 will go are determined 
uniquely. The result of Procedure 3 is an object whose second-last column is 
duplex (the lower and the upper components being c\ and C5, respectively), and 
the last column is simplex and made up of three parts: C2, C3 and C4. Thus, the 
objects produced by Procedure 3 have the following generating function: 



= (1 J 5 4%)3 -E^i("-i)-(0^) 

q t w 

- (l-qya-qt) 



Consequently, 



*w = ( i-^-# ' ft - A '"' 



q 5 t 3 w „ q 5 t 3 w . q A t 2 w ., . 



(1-^)2(1-^)3 (1-^)2(1-^)4 (l- g )2(l-^)4 

as claimed. 

We have also promised to prove the formula for A K (t). So, consider this 
procedure: 
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Figure 4: An object produced by Procedure 3. For i = 1 to 5, the cells of the 
column a arc labelled i. 



13 



Procedure 4 Given a polyomino P G S; nonempty columns c\, c 2 and C3, 
none of which has a hole; a column C4, which may be empty, but may not have 
a hole. First choose two nonadjacent cells which lie in the last column of P. 
Denote the lower chosen cell by z\ and the upper chosen cell by z 2 - Then place 
the column c\ so that the top cell of C\ becomes the right neighbour of z\, and 
place the column c 2 so that the bottom cell of c 2 becomes the right neighbour of 
z 2 . Then place the columns C3 and C4 as shown in Figure 5. 

Procedure 4 produces all elements of the set S K , but also produces some ob- 
jects which do not lie in S K . Those undesired objects are precisely the elements 
of Sg, the set whose generating function has just been found. Thus, we shall 
now compute the generating function for all objects produced by Procedure 4. 
Then we shall subtract Ag(t) and the result will be A K (t). 

In Procedure 4, if the last column of P has n cells, the cells z\ and z 2 can 
be chosen in ("2 1 ) = t" -1 ^" -2 ) ways. Once z\ and z 2 are chosen, the places 
where c\ , . . . , C4 will go are determined uniquely. The result of Procedure 4 is an 
object whose second-last column is duplex (the lower and the upper components 
being c\ and c 2 , respectively), and the last column is simplex and made up of 
two parts: C3 and C4. Thus, the objects produced by Procedure 4 have the 
following generating function: 



E^i {n ~T~ 2) ■ l(t n )A(t)} • w • (^) -j^- 

_ <l 3 tw v^°° (n-l)(n-2) Un\ A(j-\ 

~ (l-qyHl-qt)' 1 2 ' \ l l A \ l > 

q 3 tw /-1 

- (l-g)2(l_ gt )2 ' <-l- 

Consequently, 



as claimed. 



q 3 tw q 5 t 3 w 

' 61 — 7- TTTi- "tt • B\ 



(1- 9)2(1 -qty (1-^)2(1-^)3 

q 5 t 3 w qH 2 w 
+ (l-qni-qty ' Al {1-qfil-qtY ' ^ 



□ 



WehaveA(t) = A a {t)+Af}{t)+. . .+A^(t), and the expressions for A a (t) , Ap(t), A^{t) 
are given in Theorem 1. Putting these things together, we get a functional equa- 
tion for A(t). It is convenient to write that functional equation as 
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Figure 5: An object produced by Procedure 4. For i = 1 to 4, the cells of the 
column a are labelled i. 
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A(t) 



1 - qt 



1 + 



2q 3 w 



q 5 t 3 w 



(l-qni-qtf 

3qH 2 w 
{l-qf{l-qty 



1-9 



(1-gt) 2 



(1-9)2(1-^)4 



Ai 



(4) 



We solved equation ([4]) by iteration, as is usually done in the upgraded 
Temperley method. The iteration ended in 



M*) = E 



i=i (1 - q) 



E 



2i-2 



i-l„i +2i-2j.2i-l„„i-l 



(-3)-V 



<=i (1 - 9) 



2i-2 



■Ci 



+ E 



(-3) 



i-l„i 2 +4iy-2i+l 



(-3)^ i2+4i t 2i+1 u; i 



1-5 



Ai — -B] 



E 



1 (i-g) 2i - nU(i-^) 



■Ax. 



From equation (J5]), we got a system of three linear equations in three un- 
knowns: A±, B\ and C\. The first equation is just the case i = 1 of (|5]) . The 
second equation is obtained by dividing ([5]) by t, differentiating with respect to 
t and then setting t — 1. The third equation is obtained by dividing ^ by t, 
differentiating twice with respect to t, and then setting t = 1. 

Once the linear system is solved, relation (fT]) tells us how to obtain the 
sought-for generating function G. The resulting formula is stated in the next 
section. 



(5) 



4 A formula for the sought-for generating func- 
tion 

Theorem 2 The area and duplex columns generating function for simplex- 
duplex polyominoes is given by 



G(q,w) = 



NUM 
DEN ' 



(6) 



16 



where 



NUM = (1 - q) 4 (a + 7 + 2af) - 2 7 e) + q 2 w(l - qf{l+ A - ak - aft 
+ (3Z + (3X - 7K - jfj, + SZ + 5X - 2eX + 2f\Z + 2a(X - 2arjk 

— 2af}jl + 2a9X - 2f3£~X + 2~fif\Z + 2jek + 2 7 ?/I - 2 7 C? - 2^61 

— 25eX + 25fjZ) + 2q 2 w{\ - q 2 ){aX - jZ), 

DEN = (l-q) 4 (1-/3 + 5 -e + fj-aC + ae + /?? - + 1< - 7 

— 51 + Si]) - 2(1 - <7) 3 (7 + aZ) - 7 e) 

— 2q 2 w(l — q) 2 (k — (3ji + 5k — ek + (j, — £A + fjk — aC,jl + aOk 
+ (iep, - (3fjfi - j39l + f36\ + 7 C/i - lOk - 5ek + 5(Z - 5(X + 5rjk) 

— 4g 2 u>(l — q) (j3X — 7 k + a(X — af\k — (3eX + (3f\Z + jek — 7 £t) 

— 2q z w{\ — q){Z + ak — afl — f3Z + 5Z — eX + f]Z — a(X + afjk 

— afjp, + aOX + f3eX ~ [3i)Z ~ jek + 7 e/i + j(Z — 7 6>t — 5eX + 5rjZ) 
—4q a w(aX — jZ), 



E 

00 

E 

00 

E 



(sy-'q 



i-\ n V+1i-1 w i-\ 



(1 - q )2i-2 . 


"m=i(i-9 fc )" 

3 )i-y 2 +4i w i 


4 2 


(1 - g)» • 

( 


nl=i(i-9 fe )" 

_ 3 )i-l g i a +2i-2 


■(1-9*) 



^(i-g) 2< - 2 -[ni=i(i-9 fc )] -(I-? 4 ) 

'- 2 + 4-El 



2i 



fe=i 



1 _ ql J ' 



- E 



(-3) 4 -y 2+4 v 



i(i-9) 2i -[ni=i(i-9 fe )] ■ (1 _ : 



fc=l 



K = 



E 



(_ 3 )i-l g i a +2i-2 u> i-l 



i(i-9) 2i - 2 -[nt= 1 i(i-9 fe )] 4 -(i-9 i ) 



2i 



2q l 



k=i 
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-2i + 2 




In the formula for (3, one of the numbers has an overline. (Having an over- 
line does not affect the value of a number.) To obtain the formula for a from 
the formula for [3, it is enough to replace the said 2 by 1 . 

In the formula for 7, there is one 3. To obtain the formula for 5 from the 
formula for 7, it is enough to replace the 3 by 4 and change (— 3) l_1 into (—3)*. 

In the formula for £ ( resp. k ), there are two ( resp. three ) 2 's. To obtain the 
formula for e from the formula for £, and also to obtain the formula for Z from 
the formula for k, it is enough to replace each of the 2's by 1 . 

In the formula for fj (resp. \), there are two (resp. three) 3 's. To obtain 
the formula for 9 from the formula for fj, and also to obtain the formula for ft 
from the formula for X, it is enough to replace each of the 3 's by A and change 
(-3)*'" 1 into (-3)\ 

By setting w = 0, from Theorem 2 we obtain a widely known result, dis- 
covered independently by Temperley [2] and Polya [13] in the early days of 
polyomino enumeration: 

Corollary 1 The area generating function for column-convex polyominoes is 
given by 



1 - 5q + 7q 2 - 4q 3 



Proof. A short calculation. 



□ 



5 A bit of asymptotic analysis 

Formula © is very complicated, but still allows us to quickly compute Tay- 
lor polynomials of any reasonable degree. We actually computed the Taylor 
polynomial of degree 320 for the function G{q, 1). It turned out that 

G(q, 1) = q + 2q 2 + 6q 3 + 19q 4 + 63q 5 + 216q 6 + 758g 7 + 2693g 8 
+ 9608g 9 + 34269<? 10 + 121946<? n + 432701g 12 + . . . . 
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Then, for n = 2, 3, . . 
coefficient (q n ~ 1 )G(q, 1). 
91, 92, ... , 320 we have 



. , 320, we divided the coefficient (q n )G(q,l) by the 
This quotient gradually stabilizes, so that for n = 



(q n )G(g,l) 
(<r-i)G(q,l) 



3.522019812882 



rounded to 12 decimal places 



Next, for n = 1, 2, . . . , 320, we divided (q n )G{q, 1) by 3.5220198128815848301™. 
Once again, the quotient gradually stabilizes: for n — 78, 79, . . . , 320 we have 
3 5220i98"i288i584830i" = 0.119442870405. So, there is plenty of evidence that the 



coefficient (q n )G(q, 1) (i.e., the number of n-celled simplex-duplex polyominoes) 
has the asymptotic behaviour (q n )G(q, 1) - 0.119442870405-3.522019812882™. 
The dominant singularity of G(q, 1) is a simple pole, located at 0.283927988236. 

For comparison, (q n )G(q, 0) (i.e., the number of n-celled column-convex 
polyominoes) has the asymptotic behaviour (q n )G(q, 1) ~ 0.180915501882 • 
3.205569430401™. The dominant singularity of G(q, 0) is a simple pole, located 
at 0.311957055279. 



6 Further work 

Our next goal is to find the area generating function for simplex-duplex poly- 
ominoes with hexagonal cells. That should not be difficult because we now 
have a method which works with square cells. Mutatis mutandis, such methods 
usually also work when cells are hexagons. 

We are sceptical about the feasibility of an exact perimeter enumeration of 
simplex-duplex polyominoes. Anyway, it would be worthwhile to effectively do 
the computations and thus advance from scepticism to certainty. Further, it 
might be possible (although we are sceptical again) to compute the area gen- 
erating function of simplex-multiplex polyominoes. In a simplex-multiplex poly- 
omino, the columns may have any number of connected components, but it is 
not allowed that two adjacent columns have two or more connected components 
each. 

On the other hand, it is probably possible to compute the area generating 
function of simplex- duplex 2 polyominoes. Here, by a simplex-duplex2 polyomino 
we mean a column-duplex polyomino in which runs of two consecutive two- 
component columns are allowed, but it is not allowed that three consecutive 
columns have two connected components each. We feel, however, that the game 
is not worth the candle. Namely, simplex-duplex polyominoes already have a 
very complicated area generating function, and simplex-duplex2 polyominoes 
will definitely have a still more complicated area generating function. 
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